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We are interested in the phonon response in the frustrated magnets SrCr9xGai2-9xOi9 (SCGO) 
and Ba2Sn2ZnCr7xGaio-7x022 (BSZCGO). The motivation of the study is the recently discovered, 
phonon-driven, magnetic relaxation in the SCGO compound [Mutka et al. PRL 97 047203 (2006)] 
pointing out the importance of a low-energy {huo ^7 meV) phonon mode. In neutron scattering 
experiments on these compounds the phonon signal is partly masked by the magnetic signal from the 
Cr moments and we have therefore examined in detail the non-magnetic isostructural counterparts 
SrGai20i9 (SCO) and Ba2Sn2ZnGaio022 (BSZGO). Our ab-initio lattice dynamics calculations 
on SCO reveal a peak in the vibrational density of states matching with the neutron observations 
on SCO and SCGO. A strong contribution in the vibrational density of states comes from the 
partial contribution of the Ga atoms on the 2b and 12k sites, involving modes at the M-point of 
the hexagonal system. These modes comprise dynamics of the kagome planes of the pyrochlore 
slab magnetic sub-lattice, 12k sites, and therefore can drive magnetic relaxation via spin-phonon 
coupling. Both BSZCGO and BSZGO show a similar low-energy Raman peak but no corresponding 
peak in the neutron determined density of states of BSZGO is seen. However, a strong non-Debye 
enhancement of low-energy phonon response is observed. We attribute this particular feature to the 
Zn/Ga disorder on the 2d -site, already evoked earlier to affect the magnetic properties of BSZCGO. 
We propose that this disorder-induced phonon response explains the absence of a characteristic 
energy scale and the much faster magnetic relaxation observed in BSZCGO. 

PACS numbers: 74.25.Kc,78.70.Nx,78.30.-j,71.15.Mb 



I. INTRODUCTION 



The interplay between structural, electronic and dy- 
namic degrees of freedom in geometrically frustrated 
magnetic materials has the consequence of creating 
highly degenerate ground states, which have generated 
considerable interest^^^^^. Short-range correlations in 
frustrated magnets lead to the formation of weakly cou- 
pled, fluctuating clusters and consequently a macro- 
scopic, collective degeneracy can prevail^^^^^^^. 
At low temperature, frustrated magnets are expected to 
be sensitive to weak perturbations, which can raise the 
ground-state degeneracy, produce a hierarchy of closely 
space energy levels and allow the possibility of slow dy- 
namics at low temperature between the corresponding 
states. When the associated energy scales of these states 
fall well below that of magnetic interactions^, various 
perturbations including static^ or dynamic^' lattice ef- 
fects may play a key role due to coupling between lattice 
and magnetic degrees of freedom. A recent neutron spin- 
echo (NSE) examination of the pyrochlore slab antiferro- 
magnet (kagome bilayer) SrCr9a^Gai2-9icOi9 (SCGO)ii 
suggested that phonons affect the slow relaxation of this 
highly frustrated spin system. The SCGO compound 
has a particular ground-state in which dimer correlations 
within the pyrochlore slabs show a partial freezing of 
about 1/3 of the total possible ordered magnetic moment 
in spite of the high value of the average intraslab superex- 
change interaction (~ 13 meV)i^. The frozen ground- 



state does not fit into a spin-glass picture, for example 
it has been shown that both the freezing transition tem- 
perature and the frozen moment decrease with increased 
level of magnetic dilution due the Ga/Cr substitution^^. 
The extremely broad relaxation seen using NSE spec- 
troscopy was fitted with a phenomenological stretched 
exponential time-dependence to infer an activated tem- 
perature dependence with an activation energy of 
Ea ^ 7 meV^^. A vibrational mode was observed close to 
this same position by Raman spectroscopy and inelastic 
neutron scattering (INS), leading to the conclusion that 
phonons drive the relaxation. However, why phonons at 
7 meV, and not other frequencies, drive magnetic relax- 
ation is not clear and further investigation of phonons 
and related spin-phonon mechanisms is required. In this 
context we have studied the isostructural non-magnetic 
material SrGai20i9 (SGO) using ab-initio lattice dynam- 
ics calculations that provide a full picture of the phonon 
dispersion relations, total and partial density of states, 
energies of Raman active modes and the neutron scatter- 
ing cross-section. The choice of the non-magnetic coun- 
terpart SGO is motivated by experimental considerations 
(see below) and the fact that, practically, handling the 
chemical disorder associated with the Cr/Ga substitution 
in SCGO efficiently is not feasible in the lattice dynam- 
ics calculation. The calculated neutron scattering cross- 
section is used for the evaluation of the powder aver- 
aged Q-dependent intensity that can be compared with 
experiment. The calculations are accompanied by new 
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TABLE I: Lattice parameters and the relevant Ga/Cr position parameters from ab-initio calculations on SGO compared with 
experimental data on SGO and SCGO. 



parameter 


SGO (calc) 




SGO (exp)^^ 




SCGO (exp)i^ 






T=473 K 


T=298 K 


T=16 K 


T=4.2 K 


a(A) 


5.86 


5.80 


5.80 


5.79 


5.80 


c(A) 


22.87 


22.86 


22.82 


22.78 


22.66 


2a:x, y, z 


0, 0, 


0, 0, 


0, 0, 


0, 0, 


0, 0, 0^ 


2b(4e):x,2/,^^" 
12k:x,z 


0, 0, 1/4 
0.1678, 0.891 


0, 0, 1/4(0.258) 
0.1684, 0.891 


0, 0, 1/4(0.257) 
0.1683, 0.891 


0, 0, 1/4 (0.256) 
0.1682, 0.891 


0, 0, 1/4(0.256) 
0.1681, 0.892^ 




0.0278 


0.0273 


0.0273 


0.0274 


0.0283 




0.190 


0.190 


0.190 


0.190 


0.192^ 



"For refinement of experimental data it is possible to split the 2b 
occupation over two 4e sites 

^Cr occupation > 90% in least diluted samples The Kagome 
(12k) and triangular (2a) planes form the pyrochlore slab. 

^Ga occupation 100 %. 



experiments using INS and Raman spectroscopy. Exper- 
imental techniques are combined because Raman spec- 
troscopy gives a partial view of the vibrational density 
of states (F-point only) which is complemented by INS 
that can inform on the Q-dependence and characteristic 
energies throughout the Brillouin zone. 
In addition to the SCGO/SGO case we have exam- 
ined the experimental situation in the related pyrochlore 
slab compound Ba2Sn2ZnCr7a;Gaio-7x022 (BSZCGO) 
and its non-magnetic counterpart Ba2Sn2ZnGaio022 
(BSZGO). Ab-initio lattice dynamics calculations could 
not be done for BSZCGO and BSZGO due to their com- 
plex crystal structures. In BSZCGO the magnetic relax- 
ation rate at low-T was observed to be by some two orders 
of magnitude faster than in SCGO but without a well- 
defined energy scale^^. Nevertheless, both in BSZCGO 
and BSZGO Raman data reveal a phonon peak at an en- 
ergy close to the one seen in SCGO and SGO, while no 
peak is seen in the INS response. As we shall see below 
the absence of a characteristic energy in the relaxation of 
BSZCGO makes more sense when we consider the rather 
strong non-Debye enhancement of the phonon DOS that 
we attribute to the substitutional disorder, the 50/50 mix 
of Zn and Ga on the 2d site in this compound^^. Accord- 
ingly we conclude that the localized phonon modes asso- 
ciated with the substitutional disorder can also affect the 
magnetic interactions and induce relaxational dynamics 
of the magnetic system. These observations also point 
out the particularity of the frustrated magnets concern- 
ing the system dependence of the low-energy properties. 
As for SCGO, the aim of the present work is (i) to identify 
from ab-initio calculations the experimentally observed 
phonon modes (neutron and Raman spectroscopy), these 
calculations were used to obtain generalized neutron 
weighted vibrational density of states (GDOS) and the 
active Raman modes for SGO, (ii) understand the im- 
portance of non-F-point modes, and (iii) investigate the 
normal modes having an effect on the dynamics of the 
magnetic sites and hence the possibility to drive mag- 
netic relaxation. 



This paper is organized as follows: The experimental and 
computational details are provided in Section [III and Sec- 
tion [IIll respectively. Section HVl is dedicated to the pre- 
sentation of the results that are discussed in Section |Vl 
with conclusions. 



II. EXPERIMENTAL DETAILS 

Powder samples of non-magnetic SGO and BSZGO 
were prepared using a standard solid-state high tem- 
perature ceramic synthesis method and characterized 
by x-ray and neutron diffraction. As for the magnetic 
SCGO (x=0.95) and BZSCGO (x=0.97) samples, we 
analyzed those that were already used in our previ- 
ous studies^^^^^'^^^^^. The INS measurements were per- 
formed at the Institut Laue Langevin (ILL)^^ (Greno- 
ble, France), using two instruments; the cold neutron 
time-of-flight spectrometer IN6 and the thermal neutron 
time-of-flight spectrometer IN4. IN6 operating with an 
incident wavelength of = 4.1 A provides very good 
resolution (0.2 meV FWHM) in the lower energy trans- 
fer range {\huj\ < 10 meV) for the anti-Stokes spec- 
trum. On IN4 using incident wavelengths of = 2.6 
or 1.8 A one has an extended Q-range and this allows 
the Stokes spectrum to be measured at low temperature 
over a broader energy transfer range with a resolution of 
0.5 or 1 meV (FWHM), respectively. Experiments were 
performed on both SCGO/SGO and BSZCGO/BSZGO 
systems mainly at T = 300K and also at low tempera- 
tures 2 K< T <12 K for SCGO/SGO. The data analysis 
was done using ILL software tools and the GDOS was 
evaluated using standard procedures'^ without apply- 
ing multiphonon corrections. The experimental GDOS 
was normalised to 3N modes where N is the number 
of atoms in the formula unit. Raman measurements at 
T — 300 K were performed at the Institut des Materiaux 
Jean Rouxel (IMN)2£. 
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FIG. 1: (Color online) Measured S{Q,co) maps at T= 300 K for SCGO (top left) and SGO (top right). Strong paramagnetic 
scattering dominates the low-Q response of SCGO. We can note that a peaked phonon reponse comparable to the one seen in 
SGO at hcj ~ -5.5 meV with Q-dependent intensity is observed in SCGO at the highest Q-range, but at a slightly shifted energy 
transfer hu ~ -7 meV, as shown by the arrows. The bottom panels show the constant energy transfer scans ai hu — - 7 meV 
for SCGO (bottom left) and dX hix) — - 5.5 meV for SGO (bottom right). For SCGO the magnetic form factor Q-dependence 
is shown (dashed line) in order to indicate the trend of the magnetic scattering. The phonons in SCGO appear less structured 
and somewhat blurred due to the weaker coherent scattering power of Cr-atoms as well as due to the overlap with the magnetic 
signal. 



III. COMPUTATIONAL DETAILS 

Relaxed geometries and total energies were obtained 
using the projector- augmented wave (PAW) formalism^ ^ 
of the Kohn-Sham formulation of the density functional 
theory KSDFT^^^^^ at the generalized gradient approxi- 
mation level (GO A), implemented in the Vienna ab-initio 
simulation package (VASP)^^^^^. The GGA was formu- 
lated by the Perdew-Burke-Ernzerhof (PBE)^^^^^ density 
functional. All results are well converged with respect 
to A:-mesh and energy cutoff for the plane wave expan- 
sion. The break conditions for the self consistent field 
(SCF) and for the ionic relaxation loops were set to 10~^ 
eV and 10~^ eV A~^, respectively. Hellmann-Feynman 
forces following geometry optimisation were less than 
10~^ eV A~^. Full geometry optimization, including cell 
parameters, was carried out on the experimentally refined 
SGO structure^^ containing eleven crystallographically 
inequivalent atoms (5 O, 5 Ga and 1 Sr). A compar- 
ison of the ab-initio optimized and experimental struc- 
tural parameters is given in the Table [D The space 
group is PGs/mmc with 2 formula-units per unit cell 
(64 atoms). In order to determine all force constants, 
the supercell approach was used for lattice dynamics cal- 
culations. An orthorhombic supercell was constructed 
from the relaxed structure containing 8 formula-units 
(256 atoms). A second partial geometry optimization 



(fixed lattice parameters) was performed on the supercell 
in order to further minimize the residual forces. Total en- 
ergies and Hellmann-Feynman forces were calculated for 
90 structures resulting from individual displacements of 
the symmetry inequivalent atoms in the supercell, along 
the inequivalent cartesian directions (±x, ±y and ±z). 
192 phonon branches corresponding to the 64 atoms in 
the primitive cell, were extracted from subsequent cal- 
culations using the direct method^^ as implemented in 
the Phonon software^^. The coherent dynamic structure 
factor S{Q^uo) was calculated by the numerical proce- 
dure for evaluation of powder- averaged lattice dynamics, 
PALD^^'^^. For a list of wave vectors K with randomly 
chosen directions, but lengths corresponding to the Q 
range of interest, the dynamical matrix is diagonalised for 
each wave vector and the corresponding eigenfrequencies 
and eigenvectors are used to calculated S{Q^uo) for one- 
phonon creation. The wave vectors, spectral frequencies 
and intensities are then used to construct the 2-D S{Q^ uo) 
map, which can be compared with the measured map as 
will be shown below. 
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FIG. 2: (Color online) Raman spectra at T= 300 K for SGO 
(red), with calculated Raman active mode energies shown by 
the ticks at the bottom of the graph, and SCGO (blue, higher 
intensity at low energy). Both of these compounds have a 
Raman active low-energy mode in the vicinity of the kink 
seen in the GDOS, compare with Fig|3] The inset shows in 
detail the low-energy range highlighting the good match with 
the position calculated for SGO. 
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FIG. 3: (Color online) Generalized density of states (GDOS) 
at T= 300 K for SCGO (blue) and SGO (red), obtained us- 
ing the data shown over a limited ((5,a;)-range in Fig. [T] 
For SCGO the weight of the paramagnetic contribution has 
been minimized using only the high Q-data (Qei > 3.5A~^). 
Due to the remaining magnetic contribution the GDOS at 
low energy shows a faster than Debye growth as compared 
to SGO. There are distinct differences in the phonon bands 
in the range 20 meV < hco < 70 meV, but we can note a 
similar form at low energy, the peaked phonon response seen 
in SGO at hco ^ 5.5 meV and in SCGO ai hou ^ 7 meV (see 
Fig. [T]) produces a characteristic kink at that same position. 
A rather similar response of the two systems is observed in 
the high-energy range dominated by the oxygen modes. 



IV. RESULTS 



A. SGO/SCGO phonon response from experiments 

Figure [H shows the S{Q, uj) intensity maps measured at 
T = 300 K using the IN4 spectrometer, for both the mag- 
netic system SCGO and its non-magnetic counterpart 
SGO. The strong intensity at low Q centered in the elas- 
tic position in SCGO is due to the quasi-elastic magnetic 
response resulting from the fluctuating S = 3/2 moments 
of the Cr^+ atoms. Note that the Q-dependence of the 
magnetic response is modulated with a main maximum 
at Q= 1.4 and the paramagnetic form factor gives 
just an overall trend. In spite of the decay of the magnetic 
form factor with increasing Q, this signal dominates and 
swamps the phonon signal. The 5'((5, cc;)-map of SGO 
shows the phonon signal clearly visible in the absence of 
the magnetic response and it is possible to correlate the 
7 meV features at highest Q-range, Q > 3.5 A-\ m 
SCGO with the 5.5 meV features in non-magnetic sam- 
ple SGO. As we can see phonons in SCGO are masked 
due to the presence of the magnetic signal and therefore 
are difficult to measure. In addition, the chemical dis- 
order arising from the Cr substitution on the Ga sites, 
which can in principle be modelled, leads in practice to 
phonon instabilities in the calculations. For these rea- 
sons, this study is focused on the non-magnetic material 
SGO. These results are then used with proper care to ex- 



plain phonon-related phenomena in the magnetic system 
SCGO. The two systems are isostructural and accord- 
ingly the magnetic Cr sub-lattice sites can be identified 
crystallographically using the corresponding positions of 
the Ga atoms. Even though it is clear that the phonon 
response of the two systems cannot be identical both Ra- 
man spectroscopy and INS provide justification for our 
approach since the spectra of SCGO and SGO display 
similar features in the low- lying frequency range. This 
is reported in Fig. [2] for the Raman data. Both samples 
show a low-energy mode at ~ 53 cm~^ (~ 6.6 meV) and 
at - 60 cm-i (- 7.4 meV) for SGO and SCGO, respec- 
tively. Also the measured neutron density of states of the 
two systems follow quite similar trends, see Fig. [3l No 
attempt was made to eliminate the magnetic contribu- 
tion that, in spite of being small, influences the GDOS 
evaluated for SCGO in the low-energy range. Suppos- 
ing that the low-energy Raman mode are characteristic 
of Ga/Cr vibrations we can argue that the shift of po- 
sitions of the main low-frequency features is partly con- 
nected to the different masses of the Cr and Ga atoms, 52 
and 69.72, respectively. The ratio of the observed Raman 
frequencies (~0.9) is close to the square root of the ratio 
of the atomic masses (~0.86), as would be expected for 
a vibrational mode in a similar potential involving essen- 
tially the Cr or Ga atoms. Comparing the INS response 
we note a somewhat bigger difference, the peak seen at 
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FIG. 4: (Color online) Calculated and neutron determined 
generalized density of states GDOS for SCO show very sim- 
ilar overall features. The results using both IN4 (red dots) 
and IN6 (blue circles) data are in good agreement with the 
calculation (green solid line) and display a similar broadened 
peak at hu ~ 5.5 meV, see also Fig. [5] The calculated data 
has been convoluted with a gaussian to account for the ex- 
perimental resolution. 



5.5 meV in SGO appears at 7 meV in SGCO. Of course 
for modes with participation from many atomic species 
and different interatomic potentials the simple compari- 
son based on atomic mass only cannot be quantitatively 
correct. 



B. SGO/SCGO phonon response compared with 
calculation 

Fig. m compares the ab-initio determined and mea- 
sured GDOS for SGO. In order to compare with ex- 
perimental data, the calculated GDOS was determined 
as the sum of the partial vibrational densities of states 
gi{u;) weighted by the atomic scattering cross sections 
and masses: GDOS{u;)=^-{ai /Mi)gi{u;) ^ where {(Ji/Mi 
= 0.265 (O), 0.096 (Ga), 0.071 (Sr); i={0, Ga, Sr}). 
A more detailed look at the neutron data on the SGO 
compound is shown in Fig. [5] where one can see that the 
low-energy phonon response is well described by a su- 
perposition of a Debye term and a broadened peak at 
huj=5.b meV. The observed broadening of the peak {6 ^ 
2.3 meV (FWHM) does not depend on the the instrumen- 
tal resolution that is less than the apparent line-width 
both on IN4 and on IN6. Disregarding the broadening 
the calculated spectrum is in very good agreement with 
the observed one, in particular the peak at ~ 5.5 meV is 
very well reproduced and also the relative weight of the 
peak is in reasonable agreement with calculation. The 
width of the peak is not resolution limited and we con- 
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FIG. 5: (Color online) The low-energy part of the SGO- 
GDOS can be fitted (full line) with a superposition of a Debye 
term (oc u;^, Hud = 46 meV) (dashed line) and single mode 
peak with a gaussian lineshape. The results using both IN4 
and IN6 (not shown) data are in good agreement and display 
a similar broadened peak at hou ^ 5.5 meV, 6 = 2.3 meV 
(FWHM), indicating that this broadening is not due to in- 
strumental resolution. 



elude that the mode seen in Raman spectrum reported 
in Fig.[2]is included in the peak of the GDOS. Due to the 
magnetic contribution present in SCGO fitting with the 
the expected ex cj^ and a gaussian peak is not applicable. 
However, mimicking the magnetic contribution with an 
additional linear term, a reasonable fit can be obtained 
for the gaussian peak at hco = 7 meV with a width of 5 ^ 
2.4 meV (FWHM) similar to the SGO case (not shown). 
Based on the above we can consider that the experimen- 
tal data validates the computational approach in terms 
of the generalized density of states and Raman response. 
Accordingly we can go further in analyzing the calcu- 
lated partial GDOS. Fig. [6] shows total (GDOS) and the 
partials gi{u;) of SGO. The major contribution of the Ga 
atoms to the peaked low-energy response is clear. The 
Ga partial gcai^) can be further resolved into the contri- 
butions from the different crystallographic sites gcaji^)^ 
see Fig. [71 This is of interest for the case of SCGO since 
the magnetic Cr atoms occupy almost fully the sites 2a, 
12k and 4f^^ that contain only Ga in SGO, see table [H It 
appears that the 2b and 12k site partials, gGa,2b{^) and 
gGa,i2k{^) give the strongest contribution to the total 
density of states. The 2b site is fully occupied by the non- 
magnetic Ga also in SCGO while the Cr occupation on 
the 12k and 2a sites, forming the pyrochlore slab, is be- 
yond 90 % in the least diluted samples^^ (see Table |I] and 
Figure [TQ|). Structural refinements of the isostructural 
magnetoplumbite compounds, including SGO^^, usually 
indicate large Debye- Waller factor for the cation located 
in the trigonal bipyramid on the 2b site, and a split oc- 
cupation on the 4e site can be considered (see Table H]). 



6 




10 20 30 40 50 60 70 80 90 100 
Energy (meV) 



0.1 



0.05 



0.4 



0.2 




Cr/Ga 



12k 
2a 
4f . 



Ga 



2b 
4f. 



10 20 30 40 50 60 70 80 90 100 
Energy (meV) 



FIG. 6: (Color online) Calculated total and partial neutron 
weighted GDOS for SGO. Arrows indicate the peak in the 
Ga partial corresponding to the phonon mode of interest at 
hu — 5.5 meV. 



FIG. 7: (Color online) Ah initio determined partial DOS of 
Ga-atoms, gcaj, in SGO. The strongest contributions to the 
low-energy peak are from the 2b and 12k sites, respectively, 
as indicated by arrows. See Table |T] for details of site occupa- 
tions. 



The high spectral weight of the gGa,2b{^) partial is con- 
sistent with dynamical disorder being the reason for the 
large and anisotropic thermal displacement of Ga atom. 
The 12k and 2a sites constitute the pyrochlore slab mag- 
netic sub-lattice, kagome and triangular planes respec- 
tively, in the SCGO compound, and we begin to under- 
stand how phonons at 7 meV can drive magnetic relax- 
ation in SCGO. 

Calculated phonon dispersion relationships are shown in 
Fig. [8] in which the color-coded intensities correspond to 
the coherent, dynamic structure factor for one-phonon 
creation. Close to the GDOS low-energy peak (~ 6 
meV) the dispersion curves consists of flat branches with 
two modes around the M-point having maximal intensity. 
These modes make the strongest contribution to the DOS 
and we claim we see their signature in form of the gaus- 
sian peak seen in the experimental GDOS, too. At the 
F-point there are three Raman-active modes close to the 
GDOS peak but with weaker response in the INS. Due 
to the fact that the nuclei of the SGO are coherent neu- 
tron scatterers the inelastic scattering cross-section has 
characteristic features even in powder averaged phonon 
response that can also be resolved in terms of its wave- 
vector dependence. In order to gain further confidence 
in the dispersion relations and the overall phonon re- 
sponse, we have calculated the coherent structure factor 
over the whole Brillouin zone. From this data we have 
generated the S{Q^uo) map for a powder sample, using 
the recently developed, powder-averaged lattice dynam- 
ics approach^^'^^. The calculated result can be compared 
with experimental data from IN4, see Fig. [9l 

Comparison by eye of the low-energy spectrum is ham- 
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FIG. 8: (Color online) Calculated Dispersion relations along 
the high symmetry directions of the Brillouin zone in SGO, 
with color coded neutron scattering intensity, blue color cor- 
responds to the maximum intensity. The Bradley- Cracknell 
notation is used for the high-symmetry points: r= (0,0,0), 
A=(0,0,O,and M=(e,e,0). 



pered by the elastic intensity in the experimental data, 
which is not calculated. Measured and calculated S{Q^uo) 
maps are, however, in reasonably good agreement as 
shown by constant energy and constant Q cuts. 
One can also note that the modulations superposed on 
the overall Q^-dependence of the inelastic signal at en- 
ergy transfer of 5.5 meV appear periodically at positions 
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FIG. 9: (Color online) Simulated S{Q,oo) map (a) for SGO calculated for the Stokes (one phonon annihilation) process at T=0, 
Fig. [T] shows the measured anti-Stokes (phonon creation) response weighted by the thermal population. Therefore the relative 
intensity at higher absolute value of energy transfer appears less in Fig. [1] Constant Q cut (b) of measured and calculated 
S{Q,uj) (solid and dotted lines, respectively). Constant energy cut (c) of measured and calculated S{Q,(jj) ( solid and dashed 
lines, respectively). The measured cuts shown here are the phonon creation response (hu > 0) at T=2 K, note similarity of 
the constant energy cut (c) with the one shown in the lower right panel of Fig. [T]at T= 300 K. 



close to the successive M-point values at Q ~ 1.9 A~^, 
3.1 A~^, 4.3 A~^. Note that the measured data shown 
in Fig. [9] is taken at low temperature, T = 2K, and 
therefore a shorter incident wavelength (A^ = 1.8 A) was 
necessary, compared to the setting used to obtain the 
data shown in Fig. [T] (A^ = 2.6 A) , in order to reach 
high enough momentum transfer in the Stokes (neutron 
energy loss, phonon creation) side of the spectrum. The 
quality of the data at low temperature is worse due to 
reduced thermal population of phonon states. Also the 
energy resolution of the instrument is comparably worse. 
A check of the temperature dependence of the phonon 
response in SCGO in the vicinity of the spin freezing 
transition did not indicate any significant change in the 
limits of experimental precision. 

In view of the importance of the two M-point modes 
(5.6 meV and 6.2 meV), we show schematically their dis- 
placement vectors in Fig. [10] for the case of the SCGO 
magnetic compound. Therein the displacements of Cr 
and Ga atoms are highlighted. The sites 12k and 2a be- 
long to the magnetically connected pyrochlore slab net- 
work, representing the triangular planes and the kagome 
planes, respectively, and the vibrational characteristics 
there should be the most important in the context of 
the magnetic relaxation. As depicted in Fig. [10] the M- 
point phonons involve large amplitude displacements of 
the Ga atoms on the 2b sites, and the associated response 
of the 2a and 12 k sites will modulate the magnetic in- 
teractions between these sites. In this context, the ob- 
served phonon-driven magnetic relaxation could be ex- 
plained as a consequence of the strong sensivity of the 
magnetic interactions to subtle variations in the mag- 
netic exchange pathways between adjacent Cr magnetic 
cations. It was shown that the AFM coupling J between 
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FIG. 10: (Color online) Schematic representations of the nor- 
mal modes in SGO around the M-point at ^ 5.6 meV (left) 
and at ^ 6.2 meV (right), with arrows indicating the cation 
displacements. The 2a, 12k and are sites occupied in 
major part by Cr in SCGO. The 2b and 4Uv are the sites 
occupied by Ga only both in SGO and in SCGO. 



nearest-neighbour Cr ions within the pyrochlore slab is 
very sensitive to the Cr-Cr distance^. The J value fol- 
lows a phenomenological law AJ/A(i=450K/A. Because 
phonon calculations could be done for SGO only, it is 
not possible however to give an accurate value for the J 
modulation associated with the M mode in SCGO. 
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C. The case of BSZGO/BSZGCO 

We first examine ttie Raman response, Fig. [TTl Here 
again we can see that both the magnetic and non- 
magnetic compounds display a low-energy Raman peak, 
ressembling the one observed in the SGO/SCGO sys- 
tems. Also as before, one has a slightly lower energy 
for the peak in the Cr containing compound, and this 
situation can be qualitatively ascribed to the difference 
in atomic mass. The experimental GDOS has a dynamic 
range comparable to that of SGO/SCGO and especially 
the high-energy part shows similar form. Fig. [121 How- 
ever there is a qualitative difference in the experimen- 
tal GDOS as compared with SGO, the one measured for 
BSZGO does not show any distinct peak in the range 
where the Raman peaks are visible. The low-energy 
phonon response is overall rather different, it does not 
show the characteristic cj^ -dependence typical of acous- 
tic phonon response associated with the Debye-regime. 
In fact we observe an almost linear dependence GDOS 
oc uo with a progressive upwards curvature that never 
reaches a quadratic law, see Fig. [131 Such a behavior is 
quite unusual. For crystalline solids the acoustic phonons 
with linear dispersion should give a quadratic density of 
states. To explain the unconventional trend we propose 
that the substitutional disorder in the form of 50/50 oc- 
cupation by Zn or Ga on one of the 2d sites creates a 
situation that affects the lattice dynamics and gives rise 
to a disorder induced dynamic response. We can argue 
that the damping of phonon modes due to the lattice 
disorder leads to the observed non-Debye response of the 
low-energy modes, not necessarily just the long wave- 
length acoustic modes but also the ones further out in 
the Brillouin zone or even at the zone boundary. 
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FIG. 11: (Color online) Raman spectra on BSZGO (red, 
lower) and BSZSCGO (blue, upper), measured at T=300 K. 
Note the low-energy mode similar to the SGO/SCGO case 
shown in Fig. [2l 



to the thermal population nph{T). 



J pk 



g{uo)dix) 



1 



(1) 



In the case of a single mode of energy uom^ we have g{uo) oc 
S{uj — uJm)] and the phonon population and the relaxation 
rate have the activated dependence 



V. DISCUSSION AND CONCLUSIONS 

In the previous section we have provided evidence of 
the possibility of particular normal modes being respon- 
sible for the magnetic relaxation in SCGO. We suggested 
that the activated behavior of the relaxation rate is due 
to the capacity of the phonons at that energy to influence 
the magnetic system. In the following we attempt to ex- 
amine more generally the role of phonons in the temper- 
ature dependence of the magnetic relaxation. We survey 
a possible minimal model in which the important issue is 
the capability of the phonon bath to provide energy for 
the magnetic system. 



A. Magnetic relaxation controlled by phonon 
population 

Let us assume that the phonon dominated relaxation 
rate is determined by the number of occupied phonon 
states, accordingly the temperature dependence is related 



(X nph{T) 



1 



(2) 



where the approximation is valid at low temperatures 
when hoOm/kBT ^ 1. This dependence satisfies the 
experimental observation concerning the magnetic relax- 
ation in SGCO. However, if any low-energy phonon can 
induce relaxation, one can expect that the low-energy 
density of states is the important quantity, and then 
the expected T-dependence of the relaxation rate will 
depend on the functional form of the density of states 
g{uj). At low temperatures the integration in eq. [1] can 
be taken from zero to infinity and after a transformation 
of variable x = hco/kBT one finds that a power law form 
g{uj) (xujP leads to a power law in T, 



TiphiT) a {kBT/h) 



x^dx 
- 1 



(3) 



where the definite integral can be evaluated analyti- 
cally^^. Another possible situation might be such that 
there is a low-energy cut-off cOc due to e.g. Q- de- 
pendence of the effective interaction. In this case, for 
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FIG. 12: (Color online) GDOS of BSZGO (blue) compared 
with that of SGO (red), normalised as explained in the text. 
For the latter note the Debye like dependence at low temper- 
ature (dashed line with /icji) = 46 meV), see also Fig. [5] 



huJc/kBT 1, the number of phonons affecting the 

relaxation rate is evaluated to be 



nph{T) (X {kBT/hf+ 



x^dx 



(4) 

We have seen above that the experimental GDOS is a 
superposition of the Debye uj'^ and a single mode energy 
contribution in the case of SGO , while in BSZGO a linear 
cj-dependence appears. In case of independent relaxation 
processes one can expect that the relaxation rates of each 
process are additive, 



(5) 



It is straightforward to calculate nph{T) in the eqs.[T]to[4l 
Note, however, that to determine a quantitative estimate 
for the relaxation rate it would be necessary to know 
the prefactor describing the efficiency of the coupling be- 
tween the spin system and the phonon bath. Experimen- 
tally we find for the magnetic relaxation in SCGO that 
r-^ =rQ^ exp{E a/ ksT) with = 7 x 10^^ s-^ This 
means that a Debye phonon controlled relaxation of the 
type = r^^T^ would dominate at low temperatures 
for values ofr^^ > 10^^ s~^, see Fig. [HI which also shows 
that both the single mode and the cut-off picture could 
match equally well the experimental data on SCGO in the 
activated regime. It is not possible to discern any power 
law dependence either for SCGO or BSZCGO suggesting 
that long wavelength acoustic phonons are not effective 
for the magnetic relaxation process. We can also argue 
that both in SCGO and in BSZCGO at the lowest tem- 
peratures a cross-over to a quantum relaxation regime 
takes place^^'^^, and we know that at higher tempera- 
tures a quasi-elastic response with a width proportional 
to temperature prevails^^. Meanwhile, it is clear that the 
limited temperature range of the observations does not 
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FIG. 13: (Color onhne) The low-energy part of the GDOS of 
BSZGO, the solid line showing the linear non-Debye depen- 
dence for hu ^ . Data with blue and red symbols from IN6 
and IN4, respectively. 



allow a full analysis of the relaxational dynamics. We 
suggest that the most effective modes are the ones in the 
vicinity of the M-point of the reciprocal lattice, indeed 
these zone boundary modes have a spatial structure on 
the scale of the in plane lattice parameter, and there- 
fore match the nearest neighbor distance between the 
localized magnetic moments. With respect to the situa- 
tion with BSZCGO, we have observed no specific activa- 
tion energy and it is not possible to see any well-defined 
power-law dependence in the temperature dependence 
of the relaxation rate. We evoked earlier (Sect. lIVCp 
the possibility of the disorder induced modification of 
the lattice dynamics being the origin of the faster relax- 
ation, due to heavily damped phonon modes. Note that 
here again the space scale is typical of the pyrochlore 
slab in-plane unit due to the half/half Zn/Ga substitu- 
tional disorder. It is conceivable that this disorder has a 
strong influence especially on the M-point zone bound- 
ary modes that, in case of an overdamped response, can 
achieve a quasi-elastic character, and give rise to the lin- 
ear GDOS as observed, as well as to the faster relaxation 
seen in BSZCGO. Anyhow it is clear that the depen- 
dence that one might expect with reference to the model 
calculation is not observed. 



B. Other aspects of the relaxation mechanism 

Our results suggest an important role of lattice vibra- 
tions in the low-temperature spin dynamics of the frus- 
trated magnets SCGO and BSZCGO. However, one can- 
not claim that spin-phonon coupling alone could be the 
origin of the complex relaxational behavior. The sys- 
tems studied are disordered due to the non-magnetic 
dilution that has been widely studied and is a known 
factor in the pyrochlore slab compounds^^^^^— . It is 
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FIG. 14: (Color online) Modeling the relaxation, see details 
in the text. The single frequency activated model of eq. [2] 
follows the experimental data on SCGO over the range 3 K< 
T < 7 K as reported earlier and is in that temperature 
range indiscernable from the model with cut-off at the same 
frequency, eq. ID The trend of the power-law models eq. [3] 
does not match the datapoints of neither SCGO (dots) nor 
BSZCGO (squares). Note that the two models have been 
plotted with arbitrarily chosen prefactors. 

even surprising that even though the phonon properties 
of the SGO/SCGO lattices appear somewhat better de- 
fined when compared to the BSZGO, it is in SCGO that 
very strongly stretched time decay occurs. A complete 
microscopic picture of the relaxation process is still to 
be constructed and one can imagine that the fine de- 
tails of the magnetic states of lowest energy are of ma- 
jor importance in such a pursuit. Pursuing the analy- 
sis might be extremely difficult for systems like the py- 
rochlore slabs but in simpler cases details of spin-lattice 
coupling have been already examined and could pave the 
way for further work^ii^. In the general context of frus- 
trated magnetism our results point out a trend that has 
been already evoked, the low-temperature and ground- 
state properties are highly sensitive to system dependent 
perturbations. Nevertheless the present work pinpoints 
the importance of lattice vibrations and their microscopic 
character in the control of the low-energy spin dynamics. 



Ba2Sn2ZnCr7:,Gaio-7x022 (BSZCGO) based on ab- 
initio lattice dynamics simulations, inelastic neutron 
scattering and Raman measurements, in order to in- 
vestigate the origin and mechanisms of the recently 
observed^^ phonon-driven magnetic relaxation in SCGO. 
Since the magnetic signal dominates the neutron scatter- 
ing response at low energy transfer, we have performed 
new experiments on the isostructural non-magnetic 
material SrGai20i9 (SCO). Moreover, the chemical 
disorder is difficult to include in the ab-initio lattice dy- 
namics calculations so these have also focussed on SCO. 
Results of the calculations for the GDOS and S{Q, uj) for 
SGO are in good agreement with the experimental data. 
Neutron and Raman experiments show the similarity of 
the phonon response in the magnetic system with respect 
to the non- magnetic counterpart. The calculated partial 
density of states for SGO indicates that the strongest 
contribution to the total density of states stems from 
the vibrations of the Ga-atoms on the 2b and 12k sites 
of the lattice, associated with flat dispersion branches 
centered at the zone boundary M-point in the /c-space. 
Thus we conclude that the phonons response in SCGO 
corresponding to that at the M-point in SGO comprises 
the partial contribution of the atoms residing on the 12k 
sites of the magnetic sub-lattice of SCGO. The most 
relevant M-point modes have displacement vectors that 
modulate the distances between Cr sites in the kagome 
layers of the pyrochlore slab and therefore can effectively 
couple with the magnetic moments. The activation 
energy for magnetic relaxation is equal to the position 
of the characteristic peak in the experimental GDOS. 
Calculated energies of two of the Raman active T -point 
modes are included in the energy range of this peak but 
they appear with less weight in the calculated GDOS. 
In comparison BSZCGO, the other pyrochlore slab an- 
tiferromagnet shows a qualitatively different low-energy 
phonon response with an enhanced non-Debye low- 
energy response. We suggest that the faster relaxation 
without a particular energy scale is a consequence of 
this circumstance, which we attribute to the particular 
Ga/Zn disorder present in this compound. Future 
computational and theoretical work is called for to 
examine the interplay of magnetic interactions and the 
microscopic mechanisms associated with specific phonon 
displacements for better understanding of the reported 
phenomena. 



VI. SUMMARY 

We have presented a detailed investigation 
of the phonon spectra in the pyrochlore slab 
compounds SrCr9a;Gai2-9a;Oi9 (SCGO) and 
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